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Abstract 

In this paper, we study the quasineutral hmit (in other words the hmit when the 
Debye length tends to zero) of Vlasov-Poisson like equations describing the behaviour 
of ions in a plasma. We consider massless electrons, with a charge density following a 
Maxwell-Boltzmann law. For cold ions, using the relative entropy method, we derive 
the classical Isothermal Euler or the (inviscid) Shallow Water systems from fluid me- 
chanics. In a second time, we study the combined quasineutral and strong magnetic 
field regime for such plasmas. 

1 Introduction 

1.1 Physical motivation 

For high enough kinetic temperatures exceeding the atomic ionization energy, atoms tend 
to decompose into electrons and ions (that is, negatively and positively charged particles): 
a plasma is a physical or chemical system M^here such a ionization has occured. Roughly 
speaking, we simply consider that plasmas are gases made of positive and negative charges. 
Unlike gases, plasmas are highly conductive. As a consequence, particles interact with 
each other by creating their own electromagnetic fields which can dramatically affect their 
behaviour. 

The plasma state is considered as the fourth state of matter. Actually it is the most 
common one in the universe : it is widely recognized that at least 95% of the matter consists 
of it ! For instance, the suns and other stars are filled with plasma, so is the interstellar 
medium and so on. Terrestrial plasmas are also quite easy to find: they appear in flames, 
lightning or in the ionosphere. For the last decades, there has been an increasing interest 
in creating artificial plasmas, for experimental or industrial purposes. For instance, neon 
lamps or plasma displays for televisions are now part of our everyday life. An extremely 
promising application of plasmas consists in the fusion energy research (by magnetic or in- 
ertial confinement). This paper specifically aims at rigorously deriving some mathematical 
models which would help to understand the physics in tokamaks, which are the boxes in 
which plasmas from magnetic confinement fusion are contained. 



1.1.1 Basic kinetic models for plasmas 

We adopt a statistical description of the plasma: we describe the behaviour of the charged 
particles by considering kinetic equations satisfied by their repartition function. That 
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means that we do not follow particles one by one by solving Newton equations but are 
rather interested in their collective behaviour. 

We present the mathematical models we are going to study in the following. In order 
to establish them, we have to make some standard approximations which we now explain. 

• Assumption 1: We assume the plasma to be collisionless. Thus, we will consider 
Vlasov-like equations without collision operators. 

• Assumption 2: The plasma is non-relativistic and the electric field E is electrostatic. 
This means that we consider electromagnetic fields that satisfy the electrostatic ap- 
proximation of the Maxwell equations: 

r rot£^ = 0, 
I div^ = ^. 

We can also consider an additional magnetic field B satisfying 

div 5 = 



and it has to be stationary in time in order to be consistent with the electrostatic 
approximation. The electromagnetic fields act on the charged particles (with charge 
q) through the Lorentz force: 

F = q{E + v AB). 

• Assumption 3: The plasma evolves in a domain without boundaries. This means 
in particular that we may restrict to periodic data in the space variable, which may 
seem unrealistic, but which is commonly done in plasma physics and mathematics. 

We discuss the physical relevancy of these assumptions in the following Remark: 

Remark 1.1. • Discussion on Al: Let A = 47rA'^o [jf^] > ^^^re A'^o is the average 
number density of particles, eo is the vacuum permittivity, T the average temperature 
of the plasma and e is the fundamental electric charge. 

The typical collision frequency is given by: 

e^logA A^o 



So the plasma can be considered as collisionless if it is diffuse and high temperature. 
Most plasmas can be considered as collisionless to a very good approximation JWj. 

Discussion on A 2: The electrostatic approximation is relevant as soon as 

CT 

denoting by c the speed of light, r the characteristic observation time and L the 
characteristic observation length. Therefore this approximation often appears as rea- 
sonable in practical situations for terrestrial plasmas. At least, it is valid for short 
observation lengths. 
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• Discussion on A3: By ignoring boundary effects, we neglect some important physics, 
such as the formation of the Dehye sheath near walls, which are boundary layers often 
surrounding plasmas confined in some material. Very few seems to he mathematically 
known about this phenomenon, starting from a Vlasov-Poisson equation. 

Within these approximations, the kinetic system reads: 

dtfi + v.V^fi + qi/rui {E + v A B) .VJi = 
dtfe + v.V^fe + qe/rrie {E + V A B) .V^fe = 

-A.,V=i-^{qJfidv + qeJfedv), 

with X G or T" = (M/Z)" ,t; G M",t G M+. We can associate to these equations the 
initial conditions: 



fi,\t=o = fi,0: fifl > 0, / fifidvdx = 1, 

fe,\t=0 = fe,0, /e,0 > 0, f fefldvdx = I. 



[1.2) 



The parameter n is the space dimension, equal to 1, 2 or 3 in the following. Quantity fi 
(resp. /e) is interpreted as the density distribution of ions (resp. electrons) : f{t,x,v)dxdv 
is interpreted as the probabihty of finding particles at time t with position x and velocity 
V. The parameter mi (resp. rUe) is the mass of one ion (resp. electron). Likewise, qi (resp. 
Qe) is the charge of one ion (resp. electron). For simplicity we will take q^ = —e and qi = e 

We now intend to reduce the two transport equations into only one. To this end, we 
can observe that the mass ratio between electrons and ions is very smalh 

me 

— « 1, 

m 

so that qualitatively the two types of particles have really different dynamical behaviour. 
Therefore we make the additional approximation for our idealized model: 

• Assumption 4: The mass ratio between ions and electrons is infinite: ^ = +oo. 

ili'e 

This remark allows to reduce System (jl.ip to only one transport equation. Depending 
on the interpretation of Assumption 4, we get two classes of models: 

Infinite mass ions (rrij = +oo) 

One can consider the point of view of electrons, from which ions are very slow, motionless 
at equilibrium: 



"-i = y fidv = 1. (1.3) 

Then, assuming there is no magnetic field, system (II. ip written in dimensionless vari- 
ables reduces to the classical Vlasov-Poisson system: 

dtf + v.VJ + E.V,f = 

^ = -^^^ (1 4) 

f\t=o = /o > 0, / fodvdx = 1. 

This system was intensively studied in the mathematical literature (for n = 3 in partic- 
ular) and the Cauchy-Problem is rather well understood. We refer to the works of Arsenev 
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[T], Horst and Hunze |20) for global weak solutions, DiPerna and Lions for global 
renormalized solutions, Pfaffelmoser |27) . Schaeffer |30) for classical solutions, Lions and 
Perthame |22| for weak solutions with high order velocity moments and Loeper |23| on the 
uniqueness problem. 



Massless electrons (mg = 0) 

Otherwise, one can consider the viewpoint of ions : electrons then move very fast and 
quasi-instantaneously reach their local thermodynamic equilibriunj^. Then their density 
rig follows the classical Maxwell-Boltzmann law (see 



= y fedv = d{x) exp (^^r^ , (1-5) 

where V denotes the electric potential, /c^ is the Boltzmann constant, Tg the average 
temperature of the electrons, d G L^(M'^) is a term due to an external potential preventing 
the particles from going to infinity (we also refer to |3] and references therein). 
More precisely, we have: 

--g(^) 

d{x) = noe ^'bTs ^ (1.6) 

where no G K is a normalizing constant and H is the external confining potential. 
The Poisson equation then reads: 

- = j fdv - dexp (^^) . (L7) 

One should notice that in this case, in general 

fdv - dexp ^^^r^^ / 0, 

meaning that global neutrality does not hold, since the total charge of electrons is not a 
prion fixeci. 

We may also consider the case when the total charge of the electrons is fixed, in which 
case the Poisson equation reads: 

/d exp ( j^Xp 
fdv . (L8) 
/R-'^exp(^^j(irE 

The existence of global weak solutions to these two systems in dimension three has 
been investigated by Bouchut [3]. We will recall some of the properties of these solutions 
in Section [21 

An approximation widely used in plasma physics consists in linearizing the exponential 
law: 

n, = .(l + jg.). (1.9) 



^Actually since me « rrii, the typical collision frequency for the electrons is much larger than for the 
ions and thus collisions for the electrons may be not negligeable: for this reason they can reach their local 
thermodynamic equilibrium. 

^This feature will prevent us from studying this system on the torus; instead we will do so on the whole 
space R^. 
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This approximation is valid from the physical point of view as long as: 



eV 



« 1 



that is as long as the electric energy is small compared to the kinetic energy. 

We will consider this law in the case of the torus T" (with n = 1 , 2 or 3) , thus we do 
not need a confining potential (and we take d = 1). 

In the following we will only focus on models with such Maxwell-Boltzmann laws. 

1.1.2 The Debye length 

We define now the Debye length as: 



where ks is the (universal) Boltzmann constant, Tq, and are respectively the average 
temperature and density of electrons (for a = e) or ions (for a = i). 

The Debye length is a fundamental parameter which is of tremendous importance in 
plasmas. It can be interpreted as the typical length below which charge separation occurs. 

In plasmas, this length may vary by many orders of magnitude (Typical values go 
from lO^^m to 10~^m). In practical situations, for terrestrial plasmas, it is always small 
compared to the other characteristic lengths under consideration, in particular the charac- 
teristic observation length, denoted by L. Actually, the condition « L is sometimes 
required in the definition itself of a plasma. 

Therefore, if we set: 




= e << 1, 



then in many regimes, it is relevant, after considering relevant dimensionless variables, to 
consider that the Poisson equation formally reads : 



The quasineutral limit precisely consists in considering the limit e — )• 0. 
1.1.3 Why quasineutral fluid limits ? 

From the numerical point of view, kinetic equations are harder to handle than fluid equa- 
tions. Indeed the main difficulty is that we have to deal with a phase space of dimension 
6 (for x,v £ W^). Actually, another outstanding problem for simulating plasmas is the 
following : there are characteristic lengths and times of completely different magnitude 
(think of the Debye length and the observation length) that make numerics really delicate. 

In this work, we particularly aim at getting simplified hydrodynamic systems after 
taking quasineutral limits. Simplified fluid models have some advantages: 

• With a fluid description, we deal with a phase space of lesser dimension. Furthermore 
after taking the limit we now handle only one characteristic time and length. For 
these reasons, numerical simulations are easier to perform. 

Of course it is well-known that the fluid approximation is not always accurate for 
simulations of plasmas, but it is nevertheless valid in some regimes that we may 
describe in the analysis. So it is important to be aware of the physical assumptions 
we make when we derive the equations. 




(1.10) 



e^A^Ve = ±{ni- Ue) ■ 
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• Macroscopic quantities, such as charge density or current density are easier to exper- 
imentally measure (by opposition, the repartition function is out of reach). So this 
is a way to check if the initial modeling is accurate or not. 

• A simplified fluid description can help us to qualitatively describe the behaviour of 
the plasma. 

The derivation of limit models is deeply linked to issues related to the research project 
of magnetic confinement fusion. For the last few years, there has been a wider interest 
in finding simplified systems to model quasineutral plasmas for devices such as tokamaks. 
Therefore a good mathematical understanding of these becomes important, as it would 
establish some theoretical basis to compare various models, like gyrokinetic, gyrofluid, 
MHD or Euler-like equations and understand their range of validity. 

In this paper, we will focus on two kinds of quasineutral problems. In both problems, 
the starting point is the Vlasov-Poisson system with massless electrons (in other words with 
electrons following a Maxwell-Boltzmann law). First we will investigate the quasineutral 
limit alone, then we will in addition consider a large magnetic field and study the behaviour 
of the plasma in this regime. 



1.2 Quasineutral limit of the Vlasov-Poisson system with massless elec- 
trons 

In the first place, we are interested in the quasineutral limit for Vlasov-Poisson systems 
with Boltzmann-Maxwell laws and without magnetic field. We will focus in particular on 
the limit e — )■ for the following dimentionless system (we refer to the Annex for details 
on the scaling): 

System (S) : Maxwell-Boltzmann law (for x £ R^,v G M^,t G IR+) 



r dtf, + v.V^f, + E,.VJ, = 



(1.11) 



-eA,V, = J f,dv - de^^ 

. fe,\t=0 = f0,e > 0, / fo,edvdx = 1. 

The method of proof we follow allows also to treat the case of variants of system (S), so 
we will mention the results we can get, without providing complete proofs, for the systems: 

• System (S') : Maxwell-Boltzmann law with fixed total charge (for x G M.^,v G 
M3,t G M+) 

f dtf, + v.V^f, + E,.V„f, = 

Ee = -v^v; 

-eA^V, = f f,dv - ^1-12) 

, fe,\t=0 = f0,e > 0, / fo,edvdx = 1. 



• System (L) : Linearized Maxwell-Boltzmann law (for x G T^-jV G M",t G 
n = 1,2,3) 



and 



dtfe + V.V^fe + E,.V.J, = 
Ee = -VaVe 

V, - eA^V, = J hdv - 1 

fe,\t=Q = f0,e > 0, / fo,edvdx = 1. 



(1.13) 
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For systems (S) and (S') we will from now on assume the boundedness properties on 
d{x) = e--^(^): 

d = e"-^ G L^nL°^(M^). (1.14) 
V^if G W''°° for any s G N. (1.15) 



For instance, this holds for H{x) = yl + |xp. 

Remark 1.2. From the mathematical viewpoint, we have to add the confining potential H 
to ensure that the local density of electrons belongs to ( 



1.2.1 Formal derivation of the isothermal Euler system from systems (S) and 
(S') 

We will prove the local in time strong convergence of the charge density and current density: 

fedv, Je := / feVdv 



to the local strong solution (p, pu) to some Euler-type system, for initial data close (in 
some sense to be made precise later) to monokinetic data that is, 

with Uf = —. 

Let us show now how we can guess what is the limit system. First, by integrating 
the Vlasov equation against 1 and v, we straightforwardly get the local conservation laws 
satisfied by the first two moments. 

dtPe+V,.J, = 0, (1.16) 

dtJ, + V^:{fv® vf.dv] = p,E,. (1.17) 



Let us directly consider monokinetic data, i.e. f^{t,x,v) = /9e(t, x)(5^=„^(^ .j,). The local 
conservations laws reduce to: 

dtp,+Va:.(.PeUe)=0. (1.18) 



dtiPeUe) + Vx ■■ {PeUe U^) = - Pe^ x^e- (1-19) 

In the case of (S) the Poisson equation reads: 

-eA.V; = y f,dv-de^'. 

Since p^ and are uniformly bounded in L'^{L\), the following convergences hold (up 
to a subsequence) in the sense of distributions: pe ^ p and ^ J. 
If we formally pass to the limit e — )• we get: 

de^ = p. (1.20) 

Consequently, we have V = log (p/d) and therefore —VxV = — + Notice that 
= —\7xH. Thus, the asymptotic equation we can expect is the following compressible 
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Euler-type model (which can be interpreted as the isothermal Euler equation with an 
external confining force): 



dtu + u.V^u = - V:rH. 



In the case of (S') the Poisson equation reads: 



fedv 



If we formally pass to the limit e — t- we get: 

de^ 

Jde^dx ~ 



(1.21) 



:i.22) 



(1.23) 



Consequently we have V = log (^ / de^ dx) and so, we get the same Euler equation 

Remark 1.3. (Physical signification of monokinetic data) 
We define: 



Je 

V 



dvdx. 



The quantity Ti^^ is nothing but the scaled temperature of the ions. 

Considering monokinetic data corresponds to the "cold ions" assumption, that is: 



f 



J 

V 

P 



dvdx = 0, 



which means that we consider that the temperature of ions is equal to 0. 

More precisely, the cold ions approximation means from the physical point of view that 

Ti « Te. 

It turns out that this approximation is highly relevant for terrestrial plasmas and widely 
used in plasma physics, especially in the study of tokamak plasmas. 

There are two main physical reasons why it is relevant to consider that the temperature 
of electrons is much higher than the temperature of ions : first of all , there exist many 
plasma sources which can heat the electrons more strongly than the ions. Second notice 
that energy transfer in a two-body collision is much more efficient if the masses are similar. 
Thus, since ions and electrons have very different masses, there is almost no transfer of 
energy from the electrons to the ions. For instance this approximation is used in order to 
derive the classical Hasegawa-Mima equation (fl9^). 



Remark 1.4. Nevertheless we observe in the isothermal Euler limit system M.21\) that the 
ions evolve as if they had the temperature of electrons ( of order 1) ! Moreover, ions seem 
to have better confinement properties than expected, since they feel the confining potential 
in the limit equation. 



Remark 1.5. For {L) the corresponding Euler-type system is the following: 



dtp + Vx-{pu) 
dtu + u.VxU = 



0, 

-VxP- 



(1.24) 
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Actually System \l-24- ) can he interpreted in ID or 2D as an inviscid Shallow Water 



system (and p is then understood as the depth of the fluid). This quite remarkable fact is one 
amongst many analogies between geophysics and plasma physics models (see for instance 
the work of Hasegawa and Mima fl9^ and the review paper flO^). For instance the concept 
of "zonal flows " is used in both fields and the mechanism responsible for their generation 
may be the same. Only the name differs: drift waves for plasma physics, Rossby waves for 
geophysics, see Cheverry, Gallagher, Paul and Saint-Raymond for a recent mathematical 
study l^. 

1.2.2 Principle of the proof : the relative entropy method 

The relative entropy method (also referred to as the modulated energy method) was first 
introduced in kinetic theory independently by Golse [4] in order to study the conver- 
gence of solutions to a scaled Boltzmann equation to solutions of incompressible Euler for 
well-prepared data and some technical assumptions (see Saint-Raymond |29| for the latest 
developments on the topic) and by Brenier [5] in order to derive incompressible Euler equa- 
tions from the quasineutral Vlasov-Poisson equation for electrons in a fixed background of 
ions. 

More precisely Brenier shows the convergence as e — )• of the first two moments 
{pe,Je) '■= (/ fedv,f vfedv) of the starting system: 

f dtf, + v.V.J, + E,.V,f, = 

-eA,V, = ff,dv-l ^^■^^> 

. fe,\t=0 = f0,e > 0, / fo,edvdx = 1. 

to the smooth solution of the limit system which is the classical incompressible Euler 
system: 

dtU + U.Va:U + VxP = (1.26) 

div u = 0. 

Brenier treated the case of well-prepared monokinetic data (i.e. cold electrons); the 
convergence was then generalized by Masmoudi |25| for ill-prepared monokinetic data. 

We mention the works [2], [5], [7], |14| . |15| . |28| which also use the relative entropy 
method in order to derive fluid equations from Vlasov-like systems. 

The principle of the method is the following. For system (S), it can be shown that the 
following functional is non-increasing: 

jc;(t) = f^\v\'^dvdx + ^ j \V^Ve\^dx + j d{x){V, - l)e^^dx. (1.27) 

We call this functional the energy of the system. 

We then consider the functional Ti^ which is built as a modulation of this energy: 

■Heit) = ^j fe\v-u\'^dvdx + ^ j \V^V,\'^dx + j {de^' log {de^' / p) -de^' +p)dx, (1.28) 

with (p, u) a smooth solution to Isothermal Euler. 

Quite surprisingly, it turns out that the last term of Tie is similar to the usual relative 
entropy for collisional (such as Boltzmann or BGK) equations. 
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What we want to prove is that this functional is in fact a Lyapunov functional. We 
will show that indeed, He satisfies some stability estimate: 

^e(t) < ^.(0) + G,{t) + C [ \\V^u\\L^ne{s)ds, (1.29) 

with Ge{t) — T-e-s^o uniformly in time. 

Then, assuming that the initial conditions is well-prepared in the sense that 

this yields that 7ie{t) -^e^o 0, thus proving the strong convergence in some sense (which 
will be made precise later on) to smooth solutions of the isothermal Euler equation, as long 
as the latter exist. The proof relies on the fine algebraic structure of the nonlinearities in 
systems (S) and (S'). One major advantage of this method is that it only requires weak 
regularity on the solutions to the initial system but allows to prove limits in a strong sense. 
Nevertheless, it requires a good understanding of the Cauchy problem for the limit system 
(in particular, we must have a notion of stability for the limit). It should be noticed that 
even if we considered very smooth solutions (say for instance with s large) to the initial 
system, we would not be able to propagate uniform bounds and thus prove compactness. 
Indeed, the only uniform controls we have are the energy bound and the conservation of 
LP norms of the number density. 

Basically this is nothing but a stability result : roughly speaking , this result tells us 
that monokinetic solutions are stable with respect to perturbations of the energy. 

Let us mention that the method used for these two systems can also apply to the 
quasineutral limit for an isothermal Euler-Poisson version of system (S) studied by Cordier 
and Grenier [9]. We refer to Section [4.31 

1.3 Quasineutral limit for the Vlasov-Poisson equation with massless 
equations and with a strong magnetic field 

Next we are interested in the behaviour of the plasma if one applies an intense magnetic 
field. Such a regime is particularly relevant for plasmas encountered in magnetic confine- 
ment fusion research. Plasmas are expected to be confined inside tokamaks thanks to 
this magnetic field. One challenging mathematical problem is to rigorously prove if this 
strategy is likely to succeed or not. 

In this paper, we consider the simplest geometric case of a constant magnetic field with 
a fixed direction and a fixed (large) intensity. 

1.3.1 Scaling of the Vlasov equation 

We first introduce some notations: 

Notations. Let (ei,e2,e||) be a fixed orthonormal basis o/M^. 

• The subscript _L stands for the orthogonal projection on the plane (61,62), while the 
subscript \\ stands for the projection on ey . 

• For any vector X = {Xi,X2,X^\), we define X^ as the vector (Xy, —Xx,0) = X Ae\\. 
We consider a strong magnetic field of the form: 

B = Ben. 
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Roughly speaking, "strong" means that \B\ ~ 1/e. With this time a quasineutral ordering 
of the form 



Ad _ Q, 
L 



where a > is an arbitrary parameter, we get in the end the quasineutral system: 

dtfe + v.V.f, + (e, + ^) .V„/, = 



^"^^^ ,v (1-30) 

fe,\t=0 = f0,e, f fo,edvdx = 1. 

We refer to the Annex for details on the scaling. The range of parameters a > 1 is 
particularly relevant from the physical point of view. 

1.3.2 Comments on the expected result 

We will study the limit, once again by using the relative entropy method and show the 
convergence of the first two moments {p^,u^) (defined as before) to smooth solutions to the 
system: 

[ dtw + W||(9x||W = ^x^^H. 

We observe that there is no more dynamics in the orthogonal plane (that is, in the x± 
variable), which can be interpreted as a good confinement result. 

For the study of this limit we will have to face more technical difficulties than without 
magnetic field. Indeed, the strong magnetic field engenders time oscillations of order 0{l/e) 
on the number density. Consequently in order to show strong convergence we will have to: 

• filter out the time oscillations. 

• add some correction of order 0{e) to the limit {p,w) in order to get an approximate 
zero of the so-called acceleration operator (j5.13p . 

The striking point here is that we can study the limit for any value of a. In contrast, 
for the system describing the electrons with heavy ions: 

dtfe + v.VJ, + [e, + ^) .V,/, = 
E, = -VxV, 
-e''^/\xV, = ! hdv-l 

, fe,\t=0 = f0,e, f h,edvdx = 1. 

it seems primordial to take a = 1, so that the Debye length and the Larmor radius vanish 
at the same rate. This specific scaled system was studied by Golse and Saint-Raymond in 

m 

The heuristic underlying reason is that the Poisson equation with a Maxwell-Boltzmann 
law is in some sense more stable in the quasineutral limit than the "usual" one. Indeed 
the electric potential is in the limit explicitly a function of p, whereas in the "usual" case, 
it appears as a Lagrange multiplier or equivalently as a pressure. 

Remark 1.6. We could as well set ^ = S and let e,6 go to independently. One can 
readily check that we would get the same results. 
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Remark 1.7. We may also consider the linearized Maxwell-Boltzmann law and perform 
the same convergence analysis. 

1.4 Outline of the paper 

The following of this article is organized as follows. First in section [21 we recall some 
elements on the global weak solutions theory for systems (S), (S') and (L) and recall some 
useful a priori uniform bounds. This theory is due to Arsenev [1] and Bouchut [3]. Then, 
in the spirit of Lions and Perthame |22j , we consider global weak solutions satisfying some 
convenient local conservation laws. 

In section [3l we will focus on the quasineutral limit from the Vlasov-Poisson System 
(S) to an isothermal Euler system, using the relative entropy method (Theorem 13. ip . The 
crucial step is to show the algebraic identity (I3.19P that describes the decay of the relative 
entropy and from which we will be able to get a stability inequality. 

In sectionfH we generalize the method for systems (S') and (L) ( Theorems 14.11 and 14 . 2 p . 
by only sketching the proofs. We will show that this method can also be applied to the 
quasineutral limit of a system previously studied by Cordier and Grenier [9] . 

Finally, in section [5l we investigate the combined quasineutral and large magnetic field 
regime for system (S). The convergence result is stated in Theorem 15.11 

2 Global weak solutions and local conservation laws for the 
Vlasov-Poisson systems 

2.1 Global weak solutions theory 

Following Arsenev [1], it is straightforward to build global weak solutions to System (L), 
for any fixed e > 0. 

Theorem 2.1. Let n = 1, 2 or 3. We consider the functional: 

£e{t) = \j fM^dvdx + i y V.^dx + I y iV^Kpdx. (2.1) 

For any e > and initial data /o,e > bounded in L} n L°°(]R^"') such that <?e(0) is 
finite, there exists a global weak solution to (L) with E L^(L], ,,) n Lf^^^^ and £e{t) is 
non-increasing. 

Following Bouchut [3], we obtain the existence of global weak solutions to system (S) 
and (S'). We recall that d satisfies assumptions (II. 14111. 15]) (in particular, d € L^(M'^)). 

Theorem 2.2. For any e > and initial data /o,e > bounded in n L°°(M^) and 
satisfying /(I + |xp + \v\'^) fo^f^dxdv < oo: 

• The case of (S) Let J-e{t) be the functional defined as follows: 

J^^{t) = f,\v\'^dvdx + j d{x){Ve - l)e^'dx + \V^Ve\'^dx. (2.2) 

// J-"e(0) is finite, there exists f^ G L^(L^^,) n ^ global weak solution to (S) with 
J-e{t) non-increasing. 
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The case of (S'): Let Geit) be the functional defined as follows: 

Ge{t) =^ J fM^dvdx + J d{x) [Ve-\og (^j de^^dx^^ 7^^^^ 



+ 2 J l^xVel'^dx. 

If Q^{ff) is finite, there exists f^ G L'^{LI.^^) n „ global weak solution to (Sf with 
Geit) non-increasing. 

In addition, in both cases, we have : G Lf{Ll) and esssupj 3. < oo. In particular 
it means that e^' G ^Tx- 

The main difficulties in [3] are to get estimates for tlie electric potential in the Marcinkiewicz 
space to provide some strong compactness, and to use a relevant regularization scheme 
to preserve the energy inequality. 

We assume from now on that the initial data satisfy the uniform estimates : 

Ve > 0, /o,e > 0, (2.4) 
Ve > 0, /o,e G n L°°(]R"), uniformly in e, (2.5) 
3C>0,Ve>0, f,(0)<C (resp. J",, G,). (2.6) 

Using a very classical property for Vlasov equations with zero-divergence in v force 
fields, we get the following unifom in e estimates. 

Lemma 2.1. For f^ global weak solution of (L) (resp. (S), resp. (S') we have 

• (Conservation of norms) For any p £ [l,-|-oo], for any t > 0, ||/e(i)||iP < 

II/.(o)||li„- 

• (Maximum principle) If /e(0) > then for any t >0, fe{t) > 0. 

• (Bound on the energy) Vi > 0, Seit) < C (resp. J^^, resp. Ge)- 

Lemma 2.2. Define J^{t,x) = f f^vdv. Then G L'^{L]f) uniformly with respect to e. 

Proof. Actually by the same method, we can also prove that G L^{Lx) for some p > 1 
depending on the space dimension, but this result is sufficient for our purpose. The proof 
is very classical. We can first notice that there exists C > independent of e, such that: 

fe\v\^dvdx < C. 

For (L) this is clear by conservation of the energy since all the terms are non-negative. 
In the case of (S) we observe 

J fe\v\'^dvdx < F,{t) -2 j d{V, - l)e^'dx 
< J-,(0) + 2||d||ii, 

since for any a; G M, (x — l)e^ > —1. The case of (S') is of course similar. Then we can 
simply write by positivity of f^: 

< [ fe\v\dv< f fdv+ f fe\v\^dv, 



feVdv 

SO that: || Je||L^{Li) < 1 + C. 



\v\<l -'k|>l 

□ 
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2.2 Local conservation laws 



We now assume the following additional non-uniform estimate on the initial data: 

3m > 3,Ve > \v\^ h^dvdx < oo. (2.7) 

This allows to build global weak solutions with high order moments in the spirit of Lions 
and Perthame |22| . We are then able to deal with solutions satisfying local conservations 
laws for charge and current. 

For system (L), whose Poisson equation only differs from the usual one by the linear 
term V , we can straightforwardly use the results of |22) . We obtain global weak solutions 
with high order moments in v, known to satisfy as well the following local conservation 
laws: 

Lemma 2.3. Let e > 0. Let be a global weak solution to (L) with initial data satisfying 
the assumptions ^2.4^2. 7|) . Denote by p^{t^x) := / fe{t, x,v)dv and := J f^vdv. Then 
the following conservation laws hold in the distributional sense: 

dtPe + V^.J, = 0. (2.8) 
dtJe + V^:(^jv® vf.dv^ = - ^V,{V, + 1)2 



(2.9) 



+ediv^.(V^V; ® V^K) - ^VxlV^K 



This is also the case for (S) and (S'). Indeed, the only difference with Lions-Perthame's 
equations is the nonlinear term de^^ in the Poisson equation, that we may consider as a 
source term by simply noticing that de^^ £ L^ Ci L°°(M^). More precisely we get, for any 
e > 0, the (non-uniform) regularity estimates: 



Vi > 0, y \vrf,dvdx < oo, (2.10) 
GLf?,,(LinL^). (2.11) 



with A: = 1 + f . 



Thanks to these properties, we can prove that the two first local conservation laws 
(|2.12p (conservation of charge) and (|2.13p or (|2.14|) (conservation of momentum) hold in 
the sense of distributions. 

Lemma 2.4. Let e > 0. Let f^ be a global weak solution to (S) or (S') with initial 
data satisfying the assumptions ^2-4^2. 7\) . We denote the two first moments by p^{t,x) := 
J fe{t,x,v)dv and := J fevdv for any solution f^ to (S) or (S'). The local conservation 
of charge reads: 

dtPe + Vx.Je = 0. (2.12) 

The local conservation of current reads in the case of (S): 
dtJe + V^:(^jv0 vf.dv^ = -dV,(e^O + ediv,(V,K ® V,K) - |v,|V,y,|2, (2.13) 
and in the case of (S'): 

(2.14) 
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3 From (S) to Isothermal Euler 

The isothermal Euler equations (ll.2ip are hyperbolic symmetrizable. We can perform the 
change of unknown functions {p,u) i— ?• (log^,n) that leads to the system: 

r dt log(£) + V,.n + n.V, log(9 - V,F.n = 0, 

\ 9tu + n.V,.u + V,log(£) = 0. ' 

(we recall that by definition, d = e~^). 

Therefore, using classical results on hyperbolic symmetrizable systems (|24|). we get 
the local existence of smooth solutions: 

Proposition 3.1. For any initial data pQ > 0,Uo such that pQ S L^(M^), log(^) G H'^{M.^) 
and uq E H^{M.^) for s > | + 1, there is existence and uniqueness of a local smooth solution 
p > and u to ^3. 1\) such that : 

log ^, u G Ct°([o, r*[, H'{w^)) n cl{[G, t* [, h'-^{^^)) (3.2) 

pG C^([0,T*[xM3) (3.3) 

for some T* > 0. 

Since shocks may occur for large times, we will have to restrict to local results. We now 
prove the convergence of the charge and current density to the smooth solution to ()3.ip . 
as long as the latter exist. 

We restrict to well-prepared, quasi monokinetic data. 



Theorem 3.1 (The case of (S)). Let po > 0, uq verifying the assumptions of Proposition\37J 
and p, u the corresponding strong solutions of system M.21\) . We assume that the sequence 
of initial data (/o.e) satisfies the assumptions \2.Ji^2. ?]) and: 



/o,eb — uq\ dvdx — )• 0, (3.4) 

\\V~eS/.Vo4L^ ^0, (3.5) 
j {-de^o-^ log((ie^o-Vpo) - de^"'' + po) dx 0, (3.6) 
where Vo,e is solution of the nonlinear Poisson equation: 

-eA^Fo,. = j fo,edv-de^°-^. 

Then p^ weakly-* converges to p and weakly-* converges to pu in the weak sense of 
measures. Furthermore we have the following local in time strong convergences: = Je/ Pt 
strongly converges to u in the following sense: 

— Pf^dx — >• 



in and 



^fW- ^ yfp 
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Remark 3.1 (On the class of admissible initial data satisfying assumptions (j2.4ll2.7p and 
(13. 4113. 6p ). This class is not empty : indeed it includes Maxwellians of the form 



where po,e is computed by the Poisson equation after having previously chosen Vo,e such 
that ^3. 5\) and i3. 6]) hold ( for example, we can simply take Vb,e = Vq with Vq safisfying 
de^° = po and ^3. 5\) and i3. 6\) trivially hold) and Ti^^ — s-e-^o (cold ions approximation). 

Proof. Let (p, u) verifying the regularity of (13. 2113. 3p (for the moment p and u do not a 
priori satisfy the isothermal Euler equations). 

We recall that the energy for system (S) is the following functional: 

F,{t) = f,\v\'^dvdx + ^ j \V^V,\'^dx + j d{V, - l)e^'dx. 

The first two terms correspond to the energy for the quasineutral Vlasov-Poisson limit 
studied by Brenier in [5]. Therefore we accordingly modulate this quantity by considering: 



^ j f,\v-u\^dvdx + ^ j \V.^V,\^dx. 



2 

Let us now look at m. := de^^ . We can notice that: 



d{Ve — l)e^^dx = / (rrie log (m^/d) — m^)dx. 



As mentioned in the introduction, we observe a strong analogy with the relative entropy 
in collisional kinetic equations (we refer to |29) for a reference on the topic) . So by analogy, 
we modulate this quantity and hence consider the following relative entropy: 

^e(i) = ^ j fe\v — u\^dvdx + J {rUeiog {m^/ p) — + p)dx + ^ J |Vx-Ve|^da;. (3.8) 

Later on, the well known inequality (which is a plain consequence of the inequality 
X — 1 > logx, for X > 0) will be very useful: 



{^/^-Vhfdx< j {a\og{a/b) - a + b)dx. (3.9) 
We want to show that 'He(t) satisfies the inequality: 

ne{t) < UM + G,{t) + C [ \\V^u\\L^ne{s)ds, (3.10) 

Jo 

with G^{t) —7- uniformly in time. 

We show how we can deduce this kind of estimate (the computations can be rigor- 
ously justified using the local conservation laws of Lemma l2.4p . Since the energy is non- 
increasing, we have: 

dTieit) 



dt 

with: 



< hit), (3.11) 



I^{t) := / dtfe ( 2!^!^ ~ v.uj dvdx+ j f^dt f 2!^!^ ~ ^'^ ) dvdx 



(3.12) 

+ dt {me \og{d/ p)) dx + I dtpdx. 
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Let us first focus on the first two terms of /^(t). Thanks to the Vlasov equation satisfied 
by fe and after integrating by parts, we get: 

dtfe ( Trl^ip —v.u] dvdx+ I f^dt ( -\u\'^ — v.u] dvdx 



2' ' ; ./ V2 

,np — v.u 1 dvdx 
2' 



j h{dt + v.V^ + E,.V,) Qi 
j f(:{u — v).{dt + v.'Vx)udvdx — J f^E^.udvdx 

fe{u — v).{dt + u.Vx)udvdx— j fe{u — v). {{u — v).Vxu) dvdx 



PeE^.udx. 

(3.13) 

We can use the Poisson equation to compute the last term of (|3.13p . 



P^E^.udx = / Pe^xVedx 

= j dVxe^'.udx - e j ^ x ■■ (V^K ® VxVe)udx + e j ^V^| V^V^piifix 
= - J de^' divx udx - j e^'V^d-udx + e j D{u) : (V^V; (g) V^Fjdx 
- e j ^{VxYel"^ divxudx, 

where D{u) = ^ {dxiUj + dx^Ui^. . is the symmetric part of V^ii = {dxiUj)ij- 
We now focus on the last two terms of Ie{t): 

J dt{m^log{d/p)) dx + J dtpdx = j {-de^^ / p + l)dtp + J ddte^" \og{d/ p)dx 

= J{-de^'/p+l)dtp + e J dtAxVAog{d/p)dx (3.14) 
div^ J^log{d/p)dx, 



using the Poisson equation and the local conservation of mass: 

dtPe = -diVxJe- 

We observe that: 



J fe{u-v).Vx log^dvdx 



p^u.Vx^og^dx - J Je.Vx^og^dx 

de^^u.Vx log -,dx — el AajVett.Va; log —dx 
d J d 

+ / X.V^log^dx. 



(3.15) 
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Consequently, according to (|3.13p and (j3.14p we get: 



hit) = j -de^' (dtp + Vx-u + u.Vx log p)dx + j dtpdx 

+ j fe{u — v) (j)tU + U.VxU + Va; log ^ dvdx 

+e j A^y.-u.V^ log ^dx + e j dtA^VJogid/ p)dx 

f\{u - v) {{u - v).Vxu) dvdx + e / D{u) : {V^Ve VxVe)dx 



— e 



1 



VtKI^ divr udx. 



Let us now introduce the so-called acceleration operator A: 



dt log £ + Vx-u + u.Vx log £ - VxH.u 



* 5j7x + n.V,7. + V,.log£ 



We observe that : 



J p{Vx-u + uVx log p)dx= J {pVxU + u.Vxp)dx 

= J Vx-{pu)dx = 0, 



and thus f A{u,p). dx = J dtpdx. 

Gathering the pieces together we have proved: 



I,(t)= / A{u,p). 



de^' + p 



dx 



PeU - Je 

+e J AxV.u.Vx log^dx + e j dtAxV,\og{d/ p)dx 

- j fe{u - v) {{u - v).Vxu) dvdx + e J D{u) : {VxY^ V^Fjdx 

-./^|V.K|^div.n..dx. 

We are now ready to prove that satisfies the expected stability inequality. 
It is readily seen that there exists a constant independent of e such that: 



— v) {{u — v).Vxu) dvdx 



< (7 / fe\u — v\'^\\'Vxu\\L'^dvdx, 



: J D{u) : {VxV, ® VxV,)dx 
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We consider now: 

ft 



= (-^/ ^^'^e.V^ (u.V^log^) dx-ej dsV^V,.V.Aog{d/p)dx^ ds 

= -^e j V^V^.K.V^. (M.V^.log dxds + V~^J^ J VeV^V,.dsV^log{d/p)dxds 

- V^j V^^xVeit, x).V^ log{d/p{t, x))dx + j V^V^.K(0, x).V^. log(d/p(0, x))dx. 

Thanks to the conservation of the energy, -y/eV^V^ is bounded uniformly with respect 
to e in Lf^(L^). Consequently, using Cauchy-Schwarz inequality, we get for any < t < T: 

I V^('u.V^ log p/d) 11^00(^2) + ||log(p/(i)||^i,oo(^i^), 



and so we have Ge(t) — )■ when e — )• 0, locally uniformly in time. 
Finally since 

n.it) < uM + / h{s)ds, 

Jo 

we have proved that: 

n,{t) < nm + m Mu, p). (^~pfi'^j^) dx + 

+C (^J J fe\u-vf\\V,j:u\\Lo^dvdxds + e J J ^iV^^Vel^WVxuWLo^dxds^ . (3.20) 

Now we can choose p and u to be solutions of A{p, u) = 0, with initial conditions 
(/?, u)|(t=o) = (POtUq). In other words p and u are solutions to the Isothermal Euler system 

(EH). 

Then we have: 

ne{t)<nM + Ge{t) + C [ \\V ^u\\Lo^ne{s)ds (3.21) 

Jo 

and thus, since T-Le{0) — t- and Ge(t) — )■ 0, we deduce by Gronwall inequality that 

ne{t) 0, (3.22) 

when e — 7- (uniformly with respect to time). 

By inequality (13. Oh . this means in particular that 



Vde^ y/p. 



strongly in L'^L'^. 

Because of the uniform estimates in L^(L^), /o^ (resp. J^) weakly converges (up to a 
subsequence) in the weak sense of measures to some p (resp. J). In the other hand, in the 
sense of distributions, thanks to the quasineutral Poisson equation: 

de^' - Pe^O 
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and thus, 



de^^ p. 



Therefore, by uniqueness of the hmit we deduce that p = p. 

The last step of the proof rehes on a by now classical convexity argument. We first get 
the following Cauchy-Schwarz inequality: 

\J.-pM' _ {!Uv-u)dvf^ I f^\,-u\^dv. (3.23) 



Pe f fedv 

The functional (p, J) — )• J dx is convex and lower semi-continuous with respect 

to the weak convergence of measures (see [5]). Consequently the weak convergence in the 
sense of measures p^ ^ p and ^ J leads to: 

" ^""''dx < liminf / \:ll^P^dx. (3.24) 



P <:-^0 J 

So J = pu. 

To conclude, the uniqueness of the limit allows us to say that the weak convergences 
actually hold without any extraction. 

□ 

Remark 3.2 (Rate of convergence). Assume that TieiO) < C\fe. Then the previous esti- 
mates show that locally uniformly in time: 

n,{t) < C^fe. (3.25) 

4 Generalization to other quasineutral limits 

4.1 From (S') to Isothermal Euler 

Similarly, we can prove an analogous theorem for system (S'): 

Theorem 4.1 (The case of (S')). Let pQ > 0, uq verifying the assumptions of Proposition 
\3.1\ and p, u the corresponding strong solutions of system M.21\) . We assume that the 
sequence of initial data (/e,o) satisfies the assumptions \2.4-\i2. T\) and: 



strongly in L and 



JfoAv-no\^dvd.^O, (4.1) 
V^V,Vo,e ^ 0, (4.2) 



log r ^ , ^ /Po - c, V..., + Po ^ 0, (4.3) 



fde^o-.^dx ^\fde^o..^dx"^ J J de^'^^^dx 
where Vo,e is solution of the nonlinear Poisson equation: 



^0,e = j fo,i 
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Then weakly-* converges to p and weakly-* converges to pu in the weak sense 
of measures. Furthermore, we have the following local strong convergences: = Jt/ Pe 
strongly converges to pu in the following sense: 



I 



u^ — p^dx — )■ 



in and 



de^^ 



y / de^^dx 
m LT{Ll). 

Sketch of proof. According to Theorem 12. 2| the functional Qeit) is non-increasing: 
Ge{t) = \j fM^dvdx + J d{x) (^V, - log (^J de^^dx^^ J^^^^ 

+ 1 j \yxVe\^dx. 

As in the previous proof, we consider = j' de^ldx ' ^'^'^ ^® notice that: 



j d^K-log (^j de^^dx^ 



fde^'dx^^ ^ J ("^^ ('^•4) 
Since J m^dx = 1 (and thus, dt J m^dx = 0) we can actually add this quantity to the 



energy so that: 



d K - log / de^'^dx - 1 



f 

i—;-dx = / (melogfrn-e/d) — m^)dx, (4.5) 
'dx J 



fde^^ 

and we are in the same case as before. Therefore we can consider the same modulated 
energy T-Le{t). We skip the computations, which are very similar. 

□ 

4.2 From (L) to Shallow- Water 

We now treat the case of system (L), following the same methodology as before. 

4.2.1 Formal derivation of the Shallov^- Water equations 

For monokinetic data, i.e. f^{t,x,v) = Pt{t,x)6^=ue(t,x)i the conservation laws state: 

dtPe + yx-{PeU,)=Q, (4.6) 
dt{peUe) + V^. : {p,U, ® U,) = -pe^xVe. (4.7) 

We recall that the Poisson equation is: 

K - eA^y, = j f,dv - 1. (4.8) 

Since p^ and are uniformly bounded in Lf{L\), the following convergences (up to a 
subsequence) hold in the sense of distributions: p^ ^ p and J. 
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If we formally pass to the limit e — t- we get: 



V = p-l. (4.9) 

Consequently the limit system is the following: 

f dtp + Vx.{pu) = 0, ,^ 
\ dt[pu) + Vx ■ {pu 0u) = -pVxP- 

or equivalently for smooth data: 

dtp + Vx.{pu) =0, 



dfU + u.VxU = —VxP- 

As it has been said before, this system can be interpreted in ID or 2D as the Shallow 
Water equations. 

Remark 4.1 (On the kinetic version of the shallow- water limit). In a formal sense, one 
can also easily perform the kinetic limit e — ?■ and get the equation: 

dtf + v.Vxf-VxP.V„f = 0, 
p = Jfdv. 

To our knowledge, this equation is very badly mathematically understood. The only 
existence result we are able to prove is the local existence of analytic solutions. Actually 
the proof given by Mouhot and Villani J2^ (section 9, local in time interaction) in the 
Vlasov-Poisson case identically holds. Indeed, we notice that in that proof, they do not 
need the smoothing effect on the force field provided by the Poisson equation. Although it 
is not explicitly said, the case of the singular force field F = 'VxP is automatically included 
in their analysis. Of course this is not the case for the other results of their paper. 

Remark 4.2. Let us also mention that the quasineutral together with the gyrokinetic limit 
of a similar system was performed by the author in fl8^ . With quite general initial data, 
we get a limit equation of kinetic nature. In other words we do not need to restrict to 
particular initial data (or to strong regularity); this rather remarkable fact is due to the 
anisotropy of the system with the so-called finite Larmor radius scaling (}13\j). The Poisson 
equation then only degenerates in the magnetic field direction but this is overcome thanks 
to an averaging lemma. 

4.2.2 Rigorous derivation for (partially) well-prepared data 

As the Shallow Water equations (14. IIP are hyperbolic symmetrizable we get the local 
existence of smooth solutions 



Proposition 4.1. For any initial data po,Uo in i/*(T") for s > ^ + 1, there is existence 
and uniqueness of a local smooth solution to { 4.11\ ): 



p,ne C0([0, T* [, H^{T^)) n Q^([0, T* [, H^-\T^)) 

for some T* > 0. 

As before, we restrict to finite time intervals. 
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Theorem 4.2 (The case of (L)). Let po > 0,iio G iJ*(T") (s > n/2 + 1 large enough) 
and p, u the corresponding strong solutions to System We assume that the sequence 
of initial data (/o,e) satisfies the hypotheses ^2.4^2. 7\) and: 

//„,.|„-„„|^*d.^0, (413) 
/ IV^^'^Vo/'dx ^ 0, (414) 

and 

{Id - eA,)-\po,e - 1) ^ {po - 1), (4.15) 

strongly in L^. 

Then p^ weakly-* converges to p and weakly-* converges to pu in the weak sense of 
measures. Furthermore, we have the local strong convergences: 

V,^p-l 

in Lf^L"^, Ue strongly converges to u in the following sense: 

J \ ue — u\'^ Pedx — )• 0. 

Moreover, 

Remark 4.3. • Assumptions '^"-^^ ( |-<^.i5[ ) are satisfied for some smooth po as 

soon as po,e strongly converges to pQ in (For instance, when /3o,e is uniformly 
bounded in some H°' with a > 0). Indeed, in this case, we notice that: 

-^((/(i-eA,.)"Vo,.-po,e) (4.16) 

lies in a compact of H^^ endowed with its strong topology. 

This implies, thanks to the Poisson equation that y^A^^Vo^e in a compact for the 
H^^ norm. This means that -v/eV^Vo^e strongly converges to some V^^^^'o in the iP' 
norm (up to a sequence). But V^fl also strongly converges in iP' to po — 1, so by 
uniqueness of the limit in the sense of distributions, V^^'o = 0. 

• For this reason, an "ill-prepared" case would correspond to some oscillating in space 
initial data: 

Vb,e = po-l + r^, 

where r^ only weakly converges to and -v/eV^re is bounded in iP (one can think of 
r^{x) = e*^/^j. Then we would have to filter these oscillations in space to prove 
strong convergences. 

This indicates that the massless electrons have a stabilizing effect on the system, 
insofar as no time oscillations occur unless some space oscillations are imposed at 
the initial time. 

Remark 4.4. We need some additional regularity on p and u in order to handle some non- 
linear quantities : so we take s large enough (the lower bound n/2 + 1 is not sufficient). 
But we will not dwell on the optimal constant. 
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Sketch of proof . The functional £e{t) recalled below is the energy of system (L): 

The principle of the proof is to consider the following modulation of £e{t)- 

^eit) = \j fe\v - u\^dvdx + ^ J \V, - {p - dx + \V~eV^V,\^dx. (4.17) 

Then we can show with similar considerations as previous proofs that 'He(t) satisfies an 
inequality of the form: 

Tleit) <nM + Ge(.t)+ [ c\\d^u\\L^n,{s)ds 

Jo 

with Ge(t) —7- when e — ?• 0, uniformly in time. 

□ 

4.3 Quasineutral limit for the isothermal Euler-Poisson system of Cordier 
and Grenier 

As it was mentioned in the introduction, in |9], Cordier and Grenier study an isothermal 
Euler-Poisson version of (S2) and prove the quasineutral limit to the same kind of Euler 
equation (ll.2ip . So our result can be seen somehow as a generalization of theirs, since our 
startpoint is the kinetic equation. Actually the relative entropy method can also apply to 
their system. 

In [9], Cordier and Grenier consider the isothermal Euler-Poisson system (in ID): 

dtPe + dxiPeUe) = 



dtUe + UedxUe + = -OxV^ (4.18) 

—ed'i.r.Vf = Of — e 



where T is the (scaled) temperature of ions, of order 1. 

The authors perform the quasineutral limit e — )• to the so-called quasineutral Euler 
system: 

dtp + dx{pu) = 



dtu + udxu + ^^dxP = 0. 



(4.19) 



Their proof relies on rather tricky energy estimates obtained by pseudodifferential calculus 
(using the framework introduced by Grenier |17)). 

Relative entropy methods provide an alternative and more direct proof. Indeed, the 
following functional is an energy for system (j4.18p : 

<^e{t) = \j Peuldx + T j Pe(logPe " l)dx + j\v, - l)e^^ dx + \dxV,\'^dx. (4.20) 
We can consequently consider the modulated energy: 

^e{t) =^ j Pe\Ue - u\'^dx + T j Pe(log {Pe/ P) " 1 + P/Pe)dx 



+ I {melog{mJp) - me + p)dx + ^ / |5j.V;|^dx, 



(4.21) 



with = e^" . We can show, as in the previous proofs, that satisfies some stability 
inequality. In addition to the results obtained for system (S), we get the following strong 
convergence in L^(L^): 

VpT ^ Vp. (4.22) 
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5 Combined quasineutral and large magnetic field limit 



We now study the limit e — ?• of the foUowing system, which is nothing but system (S' 
with a strong magnetic field. 

dtfe + v.V^f, + (e, + ^) .V„/, = 



fe,\t=0 = f0,e, f fo,edvdx = 1. 

We shall not dwell on the existence of global weak solutions, since it is very similar to 
the theory for system (S') that was studied in section [2j 

We start with a formal analysis in order to show how we can get the expected limit 
system. 

5.1 Formal analysis 

For monokinetic data, i.e. f(:{t,x,v) = pe{t,x)S{v = u^{t,x)), the two first conservation 
laws read: 

dpe + Vx-iPeUe) = 0, 
dfUe + Ue-V xUe = E^ -\ —. 

e 

The Poisson equation reads: 

2o • 



f de^^dx 

In the limit e — t- 0, assuming that in some sense p^ ^ p^V^ ^ V (as well as V xVe 
VxV), we get: 

de^ 

^ J de^dx' ^ ^ 

and this implies that: 

p d 

If we multiply the second conservation law by e we get: 

uj- = e {dtUe + Ue-VxUe - E^) . 

This implies that 0. This convergence can not occur in a strong sense because of the 

oscillations in time of frequency 0(l/e), created by the magnetic field, but we can precisely 
describe the oscillations and consequently the strong convergence. 

We denote by TZ{s) the rotation of axis ey and angle s. Explicityly, we have: 

(cos t — sin t 0^ 
sin t cos t 
1, 

Following standard methods for singular perturbation problems (|16|. |31|). we intro- 
duce the filtered momentum defined by: 

We = n{t/e)ue. (5.3) 
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Then it is readily seen that satisfies the equation: 



dtWe + n{-t/e)w^.S/^We = TZ{t/e)E^. (5.4) 
We assume then that ^ w strongly. We take the limit e — )■ by time averaging: 

I-2TV 



1 f-^^ 



(5.5) 



and argue similarly for the other terms. We get in the end the following isothermal Euler 
system (with no dynamics in the x±_ variable): 

< VxnP (5.6) 

[ dtw + w\\dx^^w = ^x^^H. 

Of course, this system is very similar to the "usual" isothermal Euler system (ll.2ip . so 
we get the same existence result. 



Proposition 5.1. For any initial data po > 0,wo such that po £ L^{M.^), log(^) € H''' 
and wq E H^iM!^) for s > § + 1, there is existence and uniqueness of a local smooth solution 
p > and u to iS. 1]) such that : 

log ^, we C^'([0, T* [, H'{^^)) n Cl{[Q, T* [, H'-^R^)) (5.7) 
peC\[0,T*[xR-^) (5.8) 

for some T* > 0. 

From what we have seen, we can thus expect the strong convergence 

n{t/e)ue w, (5.9) 

that is to say: 

u^-TZ{-t/e)w -^0. (5.10) 

5.2 Convergence proof 

We first give the stability inequality we obtain for system (j5.ip . 



Proposition 5.2. Let (/o,e) be a sequence of initial data satisfying assumptions ^2.4^2. 7]) 
and (/e) the corresponding global weak solutions to \5.1\) . 

Lets > 3/2+1. For any sequence log ^,u^ incf{[0,T*[, H%R^))nCl{[0,T*[, H'^^iR"^)) 
we define the modulated energy: 

^e(i) = ^ / + J (l^elog^me/pe) -me+ Pe)dx + ^ j IV^V^^dx, (5.11) 

with = -j^^dx' '^^^^ following inequality holds: 

n^it) < -HeiO) + G,{t) + C [ \\VxU,\\Lo^ne{s)ds 

Jo 

+ ^7 . ).(-;!<+£<) dxd., (5.12) 
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with A^(t,x) the so-called acceleration operator defined by: 

/5tlog(|)+V..n,+n,.V.log(f)-V.i/.nA 
y dtu^ + u^.V ^u^ + V ^log [f) - -f- J 

and Ge{t) satisfying: 

G,{t) < Ce"||e°V,y,||ic«(^|) X 

(^||V^(ne-V3;logPe/d)||Lf>(L2) + \\^og{pJ d)\\^i,^ ^^j^Y (5.14) 

Proof. The proof is similar to the one given to obtain (|3.2Up in the proof of Theorem 13.11 
and therefore we omit it. □ 

As we wish to show the strong convergence (after filtering) of (p^ := J f^dv,Ue := 
f vf^dv) to solutions to system (j5.6p . a natural idea would consist in taking (pe := 
p, = TZ{—t/e)w) where p and w are the smooth solution to system (15. 6p with initial data 
Pq and wq. Unfortunately, we can not prove directly 

A,{p,n{-t/e)w) 

in a strong sense. Thus, as in |31| or |14) . we add a small correction denoted by ez^ in order 
to build a higher order approximation of the equation and make the acceleration operator 
vanish. The shape of is precisely chosen in order to "kill" the non-vanishing terms which 
only weakly converge to 0. 

In the following, we will consider the derivative of TZ: 

S{t) .= (5.15) 

which satisfies: 

( cost sint 0^ 
' - 7^(-^)_L := -sint cost 



dt 



Oy 



Theorem 5.1. Let po,W() initial data verifying the hypotheses of Proposition [3T71 with 
s > 5/2. We assume that the sequence of initial data (/o,e) satisfies the assumptions \2.4\ 
rO) and: 

ne{0)^0. (5.16) 

Let (log 2^w) the unique strong solution to ( 15. 6\) with (log ^,wq^ as initial conditions. We 
define p^ and by the relation: 

logfVf \+ev (517) 

with = [j^^ t/e)z^) ^^'^ (resp. Ze) defined by its Fourier transform J^z^ (resp. 
J^zf): 

FzP{i) = - l|^|<i^(V,^.(cS(t/e)t/;) - V ^H.S{t/e)w^) 



J^^ l|5„^|+|,|<i-^(5(t/e)^i;±)(7?).-F (v,.^ log (0) (^ - r,)dr?, 



(5.18) 
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= - l|5|<i-^ {s{t/e)V,, log 



f (5.19) 



where the operator S{t) is defined in \5.15\) . 

Then, there exists C > depending only on w and log ^, 

ll^e||L^([0,T],/i'''-i) < C, 

and locally uniformly in time we have: 

n.it) 0. (5.20) 

In particular, this means that weakly-* converges to p and weakly-* converges to 
pw^^ . Furthermore, we have the following strong convergences: 



j pe\ue-n{-t/e)w\'^dx ^ 



and 



in LfiLi 



de^^ 

Vp 



J de^^dx 



Remark 5.1. 1. Instead of a cut-off of order we could have chosen any function ^(e) 
such that for some q G]3/2, s — 1[; 



e(6 

The choice ^(e) = ^ yields a sharp convergence rate. 

2. Concerning the rate of convergence, the proof below actually shows that for any q G 
]3/2,s-l[, if: 

UM < Ce"^in(",^-9-i,i)^ 
then there exists Cq depending on q such that locally uniformly in time: 

3. When s > 7/2, we can observe in the proof that we actually do not need any cut-off 
in frequency. In this case the convergence is of order e^^^^^'°'\ 

Proof. We assume that s g]5/2,7/2]. When s > 7/2, the proof is actually much simpler, 
as we do not need any cut-off in frequency and all the estimates are straightforward. 

Step 1 We first show that ||-Ze||L^([o,T]._ff''-i) — C. Let us observe that we do not use 
the cut-off in frequency here. We have: 



/ {i + \c\'y-'\TzPAO\'d^ = m'Hs-r. 
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We then estimate: 



/ {l + \C\Y-'\T{V,,,.{S{t/e)w))\'d^ 

< Cj {l + \i\^y-'^'\FS{t/e)w\'di 

J\£.\<\ 

< c[ {l + \ey\TS{t/e)w\^d^ 

< C\\S{t/e)w\\jjs <C\\wfHs. 

Similarly we have: 

/ {l + \^\y\TiV,H.Sit/e)w^)\^dC<C\\w^\\j,.^,. (5.21) 



Finally we compute: 



<C||5(t/e)u;V,. log < C||u;||^._i || log ^||?,. 



(5.22) 



since H'^ (M"') is an algebra. This proves that there exists a constant depending on w and 
Likewise, we prove that: 

ll^f llif ([0,T],//''-i) < C". 

This yields that ez, ^ in Lf'([0, T], i^'^^). 

Now, let q g]3/2,s — 1[ be a fixed parameter. We are interested in the H'^^^ norm of 
Ze- Since q + 2 > 7/2 > s, the H'^^'^ norm of log ^ and w is not necessarily well-definecH 
and we use this time the cut-off in frequency to lower down the regularity to if ^ : 



/ [l + \i\y^'\T{V.AS{t/e)w))\''di 

J\^\<\ 

< C j {l + \i\'y{l + \i\y^'~"\FS{t/e)w\'di 

J\(.\<\ 

J\i\<^ 



Treating the other terms similarly, we obtain: 

lke||Lj°([o,T],H9+i) < C/e^'^^"''. (5.23) 

Step 2 We denote = \ ) . We introduce the filtered acceleration operator 

\w + ez^ J 

defined by: 



^When s > 7/2, q can be chosen such that q + 2 < s and thus we have ll^eHijxiQo.Ti./fi+i) ^ C*- 
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r.(- -.fdt log (^) + V^.n{-t/e)u + n{-t/e)u.V^ log (£) - V,H.n{-t/e)u\ 
^^^'""^ ^ dtu + n{-t/e)u.V,u + n{-t/e)V, log {^) J' 

(5.24) 

We show that is an approximate zero of the filtered acceleration operator in the 
sense that (we recall that q g]3/2, s — 1[): 



when e — 7- 0. 

By definition of 5, we have: 



\BeiXe)\\L^ {[0,T],H1) 0, 



dt e 

Hence we have: 

dtT{log{p/d) + ezP) 

- l\^\^iF{V.^^.{n{-t/e)w)-V.^H.n{-t/e)w^) 

- j^^ l|5„,|+|^|<i-F(7^(-^/e)7/;x(?y))^(V._, log (^))(^ - r,)di^ 

+ D,{t,i), 

where is defined by: 

A(t,0:= - el\^\<,,^(^,)F{V.,^.{S{t/e)dtw)-\I^H.S{t/e)dtw^) 

- ej^^ \_^\+\r,\<i/S{t/e)dtw^{ii)F (v,^ log (^)) (^ - n)dr^ 

- ej^^ l|5_,|+|,|<i^5(t/e)Tz;x(r?)-F(v,^9aog (^)) (i - n)dr^. 

Consequently, denoting by Bi ,, the operator in the first line of B^ (resp. i?2,e the second 
operator) we have: 



j-i?i,,(x,) = Ti,,(t,e) + T2,,(t,e) + A(t,o, 

with: 

ri,,(t,c) = i\^\^i:F{v^^.{n{-t/e)w^)-v^H.n{-t/e)w^) 

+ \.r,\+H^\Hn-t/e)w^){r,)T (v^,^ log (^)) {i - v)drj, ^^'^^^ 

T2,,{t,0 = eT {V ^.lZ{-t/e)z^ + lZ{-t/e)w.V + n{-t/e)z^ .V ^\og^ - V^H.n{-t/e)zl 
+e2^(7^(-^/e)zr.V,^f). 

(5.26) 
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Remark 5.2. Without corrector (z^ = 0) we have: 

Ti,e{t,0 = T{V^^.{n{-t/e)w^) + V^H.n{-t/e)w^) 

+ j^^ T{n{-t/e)w^){r,)T (V,^ log (^)) (C - 7?)dr? 

These terms only weakly hut not strongly converge to as e goes to ; this is why we have 
to add the corrector. 

When is defined without without cut-off in frequency, we notice that the have Ti^^{t, ^) = 

0. 

Estimating Ti 

We need the regularity of w and log p/d in order to get some decay in e for Ti^^, by 
using, for any /3 > 0: 

i,„>i<(i + ieiye2^ 



^\i-v\+\v\> 

Therefore we have: 



/ (1 + |el')'l|5|>i|-^(V.,.(7^(-^A)«;x)) M 

<(^g2.-2{g+l) f ^^^|^|2)g+l+s-(g+l)|^(^^_^yg)^^^|2^^ 
^^^2.-2{,+l)||^||2^^_ 

We handle the other terms by the same method. There exists a constant C depending 
only on q and the norm of log ^ and w such that: 

Q\i^ + \(\^y\Ti,e\^dCy <Ce'-'^-\ (5.27) 

Estimating T2^e 

We can use estimate (j5.23p . As a result, there exists a constant C depending only on 
q and the i?* norm of log ^ and w such that: 

(^j (1 + \i?y\T2,e?di^ ' < C(e X e^-'?-^ + x e^-9-2) < Ce^-^-i. (5.28) 
Estimating 

We only have dtw G ff**""*^ and consequently we do not necessarily have dtw G H'^'^^ . 
Nevertheless, we can use the cut-off in frequency to lower the regularity down to only H^~^: 

/ (1 + |eP)'^l|5|<i|-^(V.,.(5(tA)aiu;)) \'d^ 

<C [ {l + \eY+H^^^^,\T{{S{t/e)dtw))\'d^ 

^^ ,2(.+i)'-2(.-i) I, iia<i(i + \c\'r'\^mt/e)d,u^)) \'dc 

^^ ,2(g+l)-2(.-l) ll^^^llg-^- 
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Following the same method for the other terms, we finally obtain: 

1/2 



(5.29) 



Gathering the pieces together, there exists a constant C > depending on g, the 
norm of log ^, w and the norm of dt log ^, dtw such that 

As a consequence, we have proved: 



Arguing similarly for -B2,e) we finally deduce that 



(5.30) 



Step 3 Finally we check that uniformly in time : 



A,(/9e,U,). 



-me + P( 
Je - PeUe 



0. 



dxds — 0, 



as e goes to 0. We recall that and were defined in (I5.17p . 
First we have to check that 

W^eiPe, Ue)\\Lf= {[0,T],H1) 0- 

This is clear in view of Step 2, since we have: 



A(-p u)-( ^1'^^^^) ^ 



and lZ{t/e) is an isometry on any H 



We denote 



j j Al,e(/Oe,'Ue 



and evaluate: 



+ pf) dxds 



< j j \Ai^^{p^,u^)mf\dxds+ < j j \Ai^^{pf^,Ue)pe\dxds 

< C||^e(pe,?ie)||LJ«{[0,T],L°°) ( II "T-e II LJ" (LJ ) + II PII Lj°° (LJ ) II e^""" \\l° 

< C\\A,{p„u,)\\l^^[o^t],H'i) (l + ||e"'||L- 



by Sobolev embedding, H''{R'^) L°°(R^) {q > 3/2). By the estimates of Step 1, there 
exists C > independent of e such that: 



< C. 



In the other hand. 



A2,e{Pe,Ue). {Je " PeUe) dxds 

< C||A(Pe,'Ue)||L-([0,T],L-) {\\Je\\Lf(Ll) + WPeWhf (Ll)\\T^{-t h){w + eZ^)\\Lf^^ 

< Pe(Pe,«e)||Lf ([CTl.tf-J) (l + \\w\\Lf{H=) + 11^2:^11^?= 
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and the conclusions follows. 

One can also readily check, using (I5.14p . that Ge{t) — )■ uniformly in time. 
Finally this proves that 7ie{t) — )■ uniformly in time, as soon as ^e(O) — )• 0. 
Using the estimates of Step 1, we check that 

(^pe''%w + €zf^ ^ ip,w) 

in L°<=([0,T],L°°). 

In order to apply Gronwall's inequality to the inequality (I5.12p . there remains to check 
that IIViMeIIloo is uniformly bounded in e. It is sufficient to check that HcVx^^Hloo is 
uniformly bounded. According to (|5.23p and by Sobolev embedding H'^ — t- we have: 

WeV^zfh?^^ < WeV^zfh^j,. < Ce^-"-! < C. (5.31) 

Then, the other conclusions easily follow as in the end of the proof of Theorem 13. II □ 

Annex 

Scaling of the Vlasov-Poisson systems (S), (S') and (L) 

Let us introduce the dimensionless variables and unknowns: 

~ t X V 

t = — X = — V = — , 

T L Vth 

fit, X, v) = ffil i, i) V{t, x) = VV{i, x) E{t, x) = EE{i, x). 
Then the Vlasov equation with Poisson equation (II. 5p equation states: 



(5.32) 



-'-^^xVe = / fedv- edde^^' 

^ fe,\i=0 = /0,o fL^vfh I kedvdx = 1. 

In order to ensure that J fdxdv = 1, it is natural to set: 

Moreover we consider the normalizations: 

VthJ _ EL _ 
L ' y ' 



This implies that: 



Now we observe that: 



Mh = d. 

eEr _ VthT 
mvth L 



f-oV _ eokBTe _ 2 



efvf. e2 X 1/L3 
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where A^i is the Debye length. 

The quasineutral scahng consists in considering the ordering: 

with e a smah parameter. 

With this scahng we get the fohowing dimensionless system of equations (we forget the 
"for the sake of readabihty): 



. fe,\t=0 = f0,e > 0, / fo,edvdx = 1. 



(5.33) 



This system is nothing but System (S). We get Systems (S') and (L) with the same nondi- 
mensionahzation. 

Remark 5.3. To be rigorous we should also consider the confinement force —VxH on the 
ions, but we will not do so for the sake of simplicity and readability. Nevertheless we could 
handle such an external force with only minor changes in the following. 

Scaling of the Vlasov-Poisson equation (11.30ft 

We once again consider the nondimensionahzation analysis of system (S'), this time in- 
cluding the magnetic field: 



B = Be 



This yields: 



+ ^-^v.V.f, + ( + f A ei, ) .V,/, = 



EL p 



-^A^V, = efvlf f,dij-e- 



def'BTe 



(5.34) 



eV 



. fe,\i=0 = he. f^^^fh I h^dvdx = 1. 



We set = ^ : this is the cyclotron frequency (also referred to as the gyrofrequency) . 
We also consider the so-called electron Larmor radius (or electron gyroradius) r^ defined 
by: 



This quantity can be physically understood as the typical radius of the helix around axis 
ey that the particles follow, due to the intense magnetic field. 
The Vlasov equation now reads: 

d~J, + "-^nrv.Vj, + (^ftrE, + Qtv A .V J, = 0. 
L \Bvth V 

The "strong magnetic field" ordering consists in setting: 



e 



E 
Bvth 
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The quasineutral scaling we consider is: 

with a > 0. From the physical point of view, it means that we consider that both the 
Larmor radius and the Debye length vanish. We observe that: 

In most practical situations, A/j << r^, so that the range of parameters ct > 1 is 
particularly physically relevant. Finally, by having the same normalizations as before, we 
get in the end: 

dtfe + v.VJ, + [E, + ^) .V„/, = 



U\t=0 = f0,e, J fo,edvdx = 1. 
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